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ABSTRACT 


A modulo-M counter {of clock pulses) is decomposed into parallel 

modulo -m. counters, where each m. is a prime power divisor of M. Each 
1 1 a. 

m. is a cascade of a. identical modulo -p. counters, where m. = p. . The 

modulo -p^ counters are feedback shift registers which cycle through p. 

distinct states. By this organization, every possible nontrivial data frame 

subperiod (in terms of clock pulse intervals) and delayed subperiod may be 

derived. 

The number of clock pulses required to bring every (or a subset of all) 
modulo -p i counter to a respective designated state or count is determined by 
The Chinese Remainder Theorem. This corresponds to the solution of simul- 
taneous congruences over relatively prime moduli. 


vi 
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I. INTRODUCTION 


Each clock pulse interval (CPI) of a fixed-length serial data frame 
which is repetitive may be put into a one-to-one correspondence with the 
integers 0, 1, . . . , M - 1, where M is the frame length in terms of CPIs. 
Thus, a sequential network capable of assuming M distinct states can be 
used to time -tag each CPI. In effect, the sequential network is autonomous 
(i.e. , has no inputs except clock pulses) and serves as a modulo-M counter 
(of clock pulses). When synthesizing synchronous sequential logic, the 
logical designer avoids races and hazards (Ref. 1) if loading and clocking 
frequency limitations of the digital elements are respected. 

When M is not prime, every possible proper divisor (where divisors 1 
and M are excluded) corresponds to a non -trivial subperiod of the periodic 
data frame. In effect, the finite -state machine (i. e. , modulo-M counter) is 
decomposed into k submachines, where k is the number of proper divisors of 
M (Ref. 2). 

Decomposition introduces flexibility which would be costly to achieve 
in a large modulo-M counter. The combinational logic required to translate 
a state to a count is comparable for the modulo-M counter and the decom- 
posed machine. However, decomposition allows for subperiod and delayed 
subperiod generation with a modest amount of decoding logic. In a single 
machine (i.e. , modulo-M counter), the decoding logic is prohibitive for sub- 
period and delayed subperiod generation unless the count code is fixed- 
weighted. Fixed -weighted code generation in a synchronous mode, however, 
requires interstage gating (combinational logic), which grows with the 
capacity of the counter. 

Minimum overall complexity is realized with feedback shift registers 
forming submachines which operate synchronously and in series -parallel. 
FSR codes are nonweighted for every modulo (except 2). 
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II. MATHEMATICAL, BACKGROUND 


A. Fundamental Theorem of Arithmetic 

A prime number is an integer p > 1 that is divisible by 1 and p only. 

Every integer M > 1 may be uniquely expressed (except for order) as 
a prime or a product of two or more primes. The unique factorization of M 
into primes is known as the fundamental theorem of arithmetic (Ref. 3). 
That is , 


M 


“l a 2 


a 



(1) 


where the primes p^ are distinct and the exponents a are positive integers. 

B. Greatest Common Divisor 

The greatest common divisor (gcd) of two integers a and b is the 
largest integer d which divides a and b. This is denoted as 

d = ( a, b) 

where d is the largest integer, such that 

d | a and d | b 

The expression d|a (d divides a) means a = dq, where a, d and q are integers. 
Let 


a 




b 


i e 

i p 
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where a. > 0 and 3. 5 0. Then, 

l l 

min( ffj.pj) min(Qf 2 ,P 2 ) 
d = (a, b) = p x p 2 


min( a. , 3 ) 
n n 


P 


n 


Example 1 
For 


a = 588 = 2 2 • 3 • 7 2 
b = 15,435 = 3 2 • 5 • 7 3 
d = ( a, b) = 3 • 7 2 = 147 

Generally a and b are not in factored form. The Euclidean algorithm 
(Appendix I) provides an efficient means for determining (a,b) without 
employing factorization. 

A very important relationship exists between the integers a,b and 
d = (a,b), as stated in the following theorem (Ref. 3). 

If d = (a,b), there exist integers x and y, such that 

ax + by = d (2) 

An important consequence of the foregoing theorem is that if a and b 
are relatively prime, where (a,b) = 1, there exist integers x and y, such 
that 


ax + by = 1 (3) 

Conversely, if a representation such as (3) exists for 1, then (a,b) = 1. 

C . Linear Diophantine Equations and Congruences 

Diophantine equations, named in honor of the Greek mathematician 
Diophantos, are equations in one or more variables whose solutions are 
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integers (or in some cases rational numbers). Equations (2) and (3) are 
examples of linear Diophantine equations. The solution to 

ax + by = c 

where a, b, and c are given integers and x and y are integers to be deter- 
mined, involves finding an x such that ax and c yield the same remainder 
when divided by b. If a solution exists , then 

b | (c - ax) 


and y can be chosen as 


c - ax 

y = — — 

Suppose two integers s and t leave the same remainder r when divided by a 
third integer m. Then, 


s = q^ m + r 
t = q^m + r 

where 0 < r < m and s - t = (q 1 - q^) m. It follows that 

m | { s - t) 

Gauss, the 19th century German mathematician, suggested the 
following "congruence" notation: 


s = t mod m 

meaning s is congruent to t modulo m if m | (s - t). The linear Diophantine 
equation 


ax + by = c 


(4) 
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can be expressed as a linear congruence, namely 


ax = c mod b ( 5) 

Before considering the conditions for the existence of a solution or 
solutions to (4), or equivalently ( 5), let us state some additional properties 
of congruences. 

If a = c mod b and d = e mod b. Then, 

(1) a±d = c ± e mod b 

(2) ad = ce mod b 

also 

(3) ka = kc mod b for every integer k 
If ka = kc mod b and (k, b) = d, then 

(4) a = c mod (b/d) 

Note that when (k,b) = 1, then (4) becomes 

a = c mod b 

That is, the divisor k must be relatively prime to the modulo b in order to 
perform cancellation without altering the modulus. 

Returning to (4), or equivalently (5), a necessary and sufficient 
condition for a solution or solutions to exist is that 

( a, b) | c 

The number of solutions which are incongruent modulo m (hence, distinct) is 
exactly ( a, b) . 

When a and b are relatively prime, one and only one solution exists 
regardless of the integer value of c. 

Example 2 


1485x = 15 mod 2795 
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First, (1485, 2795) is determined by the Euclidean algorithm shown in 
Appendix A. 

2795 = 1485 • 1 + 1310 
1485 = 1310 • 1+175 
1310 =175-7+85 
175 = 85 • 2 + 15 
85 = 15 • 6 + 5 
15=5-3+0 

Thus, (1485, 2795) = 5 and 5 I 15. This indicates that there are five solutions 
for example 2. 

From property (4) of congruences 


1485 
5 x ~ 


-r- mod 



297x = 3 mod 559 


The foregoing congruence has one solution, since (297, 559) = 1. 
From the definition of congruences 


297x - 559y = 3 


.. _ 559y + 3 , 262y + 3 

297 = y + 297 


26 2y + 3 
297 


s (an integer) 


26 2y = 297s - 3 


y 


s + 


35s - 3 
262 
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35s - 3 

TET~ 


35s = 26 2t + 3 



— = w, u = 17w + 3 

The smallest positive value of u occurs for w = 0. It follows from w = 0, 

u = 3, t = 6, s = 45, y = 51 
and x = 96 is the solution to 

297x = 3 mod 559 

whereas 

1485x = 15 mod 2795 

has five distinct (i.e. , incongruent modulo-2795) solutions, namely, 

96, 96 + 559, 96 + 2 . 559, 96 + 3 • 559, 96 + 4 • 559 
or 96, 655, 1214, 1773, and 2332. Given the general form 

ax = c mod b, where (a,b) = d 
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A solution for x, say Xq, yields all d solutions, as follows: 


"O’ 0 


+ — 

+ d’ 


. 2b 
X 0 + d * 


x 0 + 


(d - l)b 


Corresponding to each solution for x is a solution for y in the linear 
Diophantine form 


ax - by = c 

However, as in example 2, values of y are often not required. 

A more efficient method for finding a solution for x is given in 
Appendix I. 

D. The Chinese Remainder Theorem for Integers 

The Chinese Remainder Theorem guarantees a unique solution for 
simultaneous congruences over moduli which are relatively prime by pairs. 
The theorem may be stated as follows (Ref. 3): 

Every system of linear congruences in which the moduli are relatively 
prime in pairs is solvable, the solution being unique modulo, the product of 
the moduli. 

Given the simultaneous congruences 

x = a^ mod 
x = a_ mod m_ 


x = a mod m 
n n 


( 6 ) 


8 


JPL Technical Memorandum 33-539 



where (m. , m.) = 1 for all i, j, where i j and a 1 , a,, • • • , a are any set of 
i J 1 c n 

integers, let 


and 


M = rrijin^ 


m 


n 


M. 


M 


i m. 

l 


Since (M^,rm) = 1, a unique solution exists for in the linear congruence 


M.y. = 1 mod m. for all i 
i 7 i i 

There is one and only one solution for x, which is determined as 
follows: 


n 


-Evt 

i=l 


M. 
i 


mod M 


(7) 


Note that, as expressed in (7), x is a solution of each congruence in (6). 

a.y.M. = a. mod m. 

1/1 ii l 

= 0 mod m . , where i ^ i 
J 

The latter results since m. is a factor of Ivh . The value of x is such that 
0 < x < M. 

Example 3 


x = 1 mod 3 a^ = 1 

x = 2 mod 4 a^ = 2 

x = 3 mod 5 a^ = 3 


mj =3 
m 2 = 4 

m 3 = 5 
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M = 3 • 4 • 5 = 60 


Mj = 20, M 2 = 15, M 3 = 12 
20y ^ = 1 mod 3 
1 5y ^ = 1 mod 4 
12y 3 = 1 mod 5 

Unique solutions for y ^^2 an< ^ ^3 are ^ an< ^ res P ec tively . 

x = (40a^ + 45a 2 + 36a^) mod 60 
x = (40 • 1 + 45 • 2 + 36 • 3) mod 60 
x =58 mod 60 


Check 

58 =1 mod 3 
58 =2 mod 4 
58 = 3 mod 5 

A modulo-3, a modulo-4, and a modulo-5 counter would be in state 12 3 
(i. e. , a^ = 1, a 2 = 2 and a^ = 3) for n = 58 + k60 CPI, where k = 0, 1, 2, • • • 
(State 12 3 repeats every 60 CPIs; see Table 1.) 


10 
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in. APPLICATION OF THE CHINESE REMAINDER THEOREM 
TO TIMING AND CONTROL SIGNAL GENERATION 

The following data frame length has been proposed for the Mariner 
Venus -Mercury 1973 spacecraft: 

M = 14, 817, 600 

Unique factorization (except for order) gives 

,, ,6 ~3 c 2 _3 

M = 2 3 5 7 


Let 


M 


m. 


m 2 m 3 


m 


where 

= 2^ = 64 

m 2 = 3 3 = 27 

= 5^ = 25 

m 4 = 7 3 = 343 

The arrangement of rm factors is arbitrary as long as pairwise relative 
primeness holds. The unique factorization into products of powers of 
distinct primes, where each distinct prime power corresponds to an m. , 
guarantees pairwise relative primeness. Furthermore, prime power factori- 
zation enables one to enumerate all the divisors of M. Let 


P 1 = 2 > P 2 = 3 » P 3 = 5 and P 4 = 7 
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The following polynomial factors, when multiplied, yield terms which 
correspond to every divisor of M: 


( 1 + Pj + Pi + Pj + + Pj + Pj) 

x ( l + p 2 + P2 + p|)( 1 + P 3 + P3) 

XU + P 4 + P4 + P4) 

The number of terms in the resulting polynomial is 

(6 + 1)( 3 + 1)<2 + 1)(3 + 1) = 336 

Thus, there are 336 divisors of M, 334 of which are proper (1 and M are 
improper). 

As shown in Figure 1, there are four counters, which count in modulo 
m^, modulo m^, modulo m^, and modulo m^. Each is composed of identical 
cascaded feedback shift registers (FSRs), where the number of states a 
particular FSR cycles through is a prime factor of M. Each FSR is designed 
to operate synchronously. Furthermore, the successive states through which 
a modulo-p^ FSR cycles (except for p^ = 2) is not in binary order. Unused 
states are always driven into the major cycle (see Ref. 2). An output cor- 
responding to a scale of p. is clocked to the next modulo -p. counter to realize 

2 1 1 
a modulo -p. count, etc. 

As shown in Fig. 1, all rm states associated with the modulo-im 
counter may be decoded. In practice, however, only those states corre- 
sponding to a count of interest are decoded. 

Example 4 

Assume that CPI 10,942 must be identified. 

10,942 = 62 mod 64 

= 7 mod 27 

= 17 mod 25 

= 309 mod 343 


12 
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CPI 10, 942 occurs when , a^, a^, and are 62, 7, 17, and 309, respectively. 
The occurrence is once per frame. That is. 


CPI x = 10, 942 + kM for k = 0, !,••• 


Example 5 

CPI x = 0 + k400 for k = 0, 1, • • • is to be generated. That is, M is 
to be divided into equal subframes of length 400. 


Let 


4 2 

400 =2-5 



M = m^m^ = 400 


The content of the leftmost four stages of m^ and all six stages of m^ must be 
detected, so that when all zeros are stored, an output x is generated. The 

a 

first output may be delayed up through M - 1 CPIs (up to 399) by detecting an 
appropriate nonzero a^a ^ combination. Assume that the first output is to be 
delayed 72 CPIs. 


7 2=8 mod 1 6 
= 22 mod 25 

The combination a^a^ ° f 8 22 will appear at CPI x = 72 + k400, where 


k = 0, 


1 , 
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Example 6 


In example 5, assume that a^ = 15 and a^ = 20. The number of CPIs 
designated by x required for this a^a^ combination to appear can be deter- 
mined by the Chinese Remainder Theorem. 


x =15 mod 16 
x = 20 mod 25 


M 


1 


M _ 400 oc 

m^ " HT = 25 


A 

M 


2 



16 


25y^ s 1 mod 16 
1 6y^ = 1 mod 25 


Unique solutions for y 1 and y 2 are 9 and 11, respectively. 

x = (255a^ + mod 400 

x = 95 + k400 for k = 0, 1, 


Note that 


9 5 = 15 mod 16 
== 20 mod 25 


Thus, the a^a^ combination (i.e. , count) of 15 20 will appear for the first 
time at CPI 95 and will reappear every 400 CPIs thereafter. 

Example 7 

Given 

M = 2 6 3 3 5 2 7 3 = 14,817,600 
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where 


rrij = 2^ = 64 

m 2 = 3 3 = 27 

= 5 3 = 25 

m . = 7 3 = 343 

4 

and 

M = — = 231, 525 

1 m i 

M 0 = — = 548, 800 

2 m 2 

M_ = — = 592, 704 

3 m 3 

M = = 43, 200 

m 4 

The minimum number of clock pulses required to realize the counts 

a, = 37, a_ = 17, a_ = 2, a. = 341 
1 2 3 4 

is determined as follows: 

A unique solution for each y\ is guaranteed, where 

231, 525 y = 1 mod 64 

548,800 y 2 = 1 mod 27 

592, 704 y 3 he 1 mod 25 

43, 200 y 4 = 1 mod 343 
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The solutions are 

y i = 45 > y 2 = 13 > y 3 = y 4 = 

and 

x s + M^y 2 a 2 + M 3^3 a 3 + M 4 y 4 a 4) mod M 

x = 9, 039, 077 + kl4, 817, 600 for k = 0,1, 2,--- 

The minimum number of clock pulses required is 

x = 9,039,077 clock pulses 


Check: 


9,039,077 = 37 mod 64 

= 17 mod 27 

5 2 mod 25 

= 341 mod 343 
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IV. RECURRENCE RELATIONSHIPS FOR THE MODULO-p. COUNTERS 

An r -stage feedback shift register which provides a count modulo -p^ 
may be characterized by an rth-order recurrence relation (i. e. , an rth-order 
difference equation). 


b k “ f(b k-l’ b k-2’ 


, h. ) 

k-r 


( 8 ) 


The bit fed back at CPI k is denoted by b^. The content of the ith stage at 
CPI k becomes the content of the {i + l)th stage at CPI k + 1. That is, 

b k-i = b (k+l)-(i+l) ( 

Expression (9) accounts for the shifting in the register. From (8), the bit 
being fed back is a Boolean function of the contents of the register. The 
initial state of stage i is b where CPI k = 0. The number of stages 
required for a modulo-p^ counter is r, where r satisfies the following 
inequalities : 


,r-l _r 

2 < p. < 2 


The recurrence relationships for the proposed modulo-p. counters in 
example 7 and Fig. 1, where p^ = 2 , p^ = 3, p^ = 5, and p^ = 7, are: 

P 1 = 2 b k = 1 © b k _i = b k _! 


P? = 3 V = 


b k-l b k-2 


p 3 


= 5 b, = 


b k-2 b k-3 


P 4 7 b k b k-l b k-2 b k-3 + b k-l b k-3 


( 10 ) 
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The symbols (+}, 1 , and + denote exclusive— or complementation and 
inclusive— or (logical summation), respectively. Juxtaposition denotes 
logical multiplication. The initial state of each stage is assumed to be 0 
(i. e. , b ^ = 0 for all i). 

State diagrams of each of the four FSRs appear in Fig. 2. Note that 
unused states are always driven into the desired cycle of states. 

The FSRs appearing in Fig. 1 are not shown in detail. In particular, 
the feedback networks are not explicitly drawn. See Fig. 3 for a generalized 
FSR (without a decoding network). The state of the ith stage at clock pulse 
interval (CPI) is denoted by b^ The binary digit being fed back at CPI k 
is denoted by b^. 

The content of the first or leftmost stage is replenished by b, after 

K 

shifting, where b. is a Boolean function of b. . , b. b, , . , b, . The 

k k-1 k-2 k-r+1 k-r 

right-hand side of the recurrence relations in (10) corresponds to Boolean 
feedback functions. 

FSR implementations for (10) appear in Fig. 4. The shift register 
stages are 1-enable JK flip-flops , whose characteristic equation is 


Q = Jq' + K'q 


where J and K are 1 -enable inputs and q and Q are the present and next state, 
respectively. 

Symbolism in Fig. 4 has been simplified, using the following 
correspondences: 


b. . <— * b. and b, b 
k-i i k 


J 


k-i 


J. 

l 



K. 

l 
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Table 1. A modulo-60 counter decomposed into parallel modulo m^ = 3, 
modulo mg = 4, and modulo mg = 5 counters 


CPI 

n 

Count state 

CPI 

n 

Count state 

m i 

m 2 

m^ 

m l 

m 2 

m 3 

0 

0 

0 

0 

30 

0 

2 

0 

1 

1 

1 

1 

3 1 

1 

3 

1 

2 

2 

2 

2 

32 

2 

0 

2 

3 

0 

3 

3 

33 

0 

1 

3 

4 

1 

0 

4 3 

34 

1 

2 

4 

5 

2 

1 

0 

35 

2 

3 

0 

6 

0 

2 

1 

36 

0 

0 

1 

7 

1 

3 

2 

37 

1 

1 

2 

8 

2 

0 

3 

38 

2 

2 

3 

9 

0 

1 

4 

39 

0 

3 

4 

10 

1 

2 

0 

40 

1 

0 

0 

11 

2 

3 

1 

41 

2 

1 

1 

12 

0 

0 

2 

42 

0 

2 

2 

13 

1 

1 

3 

43 

1 

3 

3 

14 

2 

2 

4 

44 

2 

0 

4 

15 

0 

3 

0 

45 

0 

1 

0 

16 

1 

0 

1 

46 

1 

2 

1 

17 

2 

1 

2 

47 

2 

3 

2 

18 

0 

2 

3 

48 

0 

0 

3 

19 

1 

3 

4 

49 

1 

1 

4 

20 

2 

0 

0 

50 

2 

2 

0 

21 

0 

1 

1 

51 

0 

3 

1 

22 

1 

2 

2 

52 

1 

0 

2 

23 

2 

3 

3 

53 

2 

1 

3 

24 

0 

0 

4 

54 

0 

2 

4 

25 

1 

1 

0 

55 

1 

3 

0 

26 

2 

2 

1 

56 

2 

0 

1 

27 

0 

3 

2 

57 

0 

1 

2 

28 

1 

0 

3 

58 

1 

2 

3 

29 

2 

1 

4 

59 

2 

3 

4 
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2 CASCADED 5-STATE 
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7-STATE 

FSRs 


Fig. 1. Proposed spacecraft timing and control signal generation for 
Mariner Venus -Mercury 1973 
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Fig. 4. Implementation of FSRs whose state diagrams appear in Fig. 2 
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APPENDIX A. ALGORITHMS FOR FINDING THE GREATEST 

COMMON DIVISOR 


I. THE EUCLIDEAN ALGORITHM 

The greatest common divisor of two integers may be computed by means 
of the Euclidean Algorithm (Refs. 3 and 4). 

Let r 2 and r ^ denote two positive integers, where 


r -2 r -l h 0 + r 0 


r -l = r 0 h l + r l 


r k r k+l h k+2 + r k+2 


0< r 0 < r -l 


0< r l < r Q 


° < r k+2 < r k+l 


r = r h + r 
n-3 n-2 n-1 n-1 


r _ = r ,h + 0 
n-2 n-1 n 


0 < 


n-1 


r 

n 


2 


0 = r < r , 
n n-1 


( 11 ) 


Note that successive values of r are decreasing integers, so that r^ = 0 for 
some integer n. Also 


^ r -2 ,r -P (r -r r o } “ ^ r o ,r i^ " 


= (r ) = (r , 0 ) = r 

n-2 n-1 n-1 n-1 


Linear Diophantine equations, and hence linear congruences, can be 
solved by means of the Euclidean Algorithm. This can best be illustrated by 
an example. 
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Example 8 


The linear congruence in example 7, 

231, 525 = 1 mod 64 

was asserted to have the solution y^ = 45. The congruence may be expressed 
as a linear Diophantine equation as follows: 

231,525 y^ - 64z = 1 

where y^ and z are integers. First, the Euclidean Algorithm is employed to 
determine (231525, 64). 

231,525 = 64 • 3617 + 37 


64 

= 37 

• 1 

+ 

27 

37 

= 27 

• 1 

+ 

10 

27 

= 10 

• 2 

+ 

7 

10 

= 7 ■ 

1 + 

3 


7 

= 3 • 

2 4 

1 


3 

= 1 ' 

3 




and (231525, 64) = 1. Note that h Q = 3617, h 1 = 1, h.^ = 1, h 3 = 2, h.^ = 1 , and 
h,- = h n _j = 2. Also, r^ = r ^ = 1 =(231525,64). Using successive partial 
results of the Euclidean Algorithm in reverse yields 

1 = 7 - 3 • 2 

= 7 - (10 - 7 • 1)2 = 7 • 3 - 10 • 2 

} 

= (27 - 10 • 2)3 - 10 • 2 = 27 • 3 - 10 • 8 
= 27 • 3 - (37 - 27 • 1)8 = 27 • 11-37-8 
= (64 - 37 • 1)11 - 37 • 8 = 64 • 11 - 37 ■ 19 
= 64 • 11 - (231, 525 - 64 • 3617)19 
= 64 • 68,734 - 231, 525 • 19 
= 231,525 (-19) + 64 • 68,734 
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Thus , 


and 


= -19 = 45 mod 64 


z = 68,734 


The foregoing method is inefficient, since all the partial results of the 
Euclidean Algorithm must be stored in order to solve and z. Further- 
more, in solving a linear congruence, the value of z (i.e. , the multiple of 
the modulo) is not required. 

II. SUCCESSIVE CONVERGENTS 

The Euclidean Algorithm, together with two additional recurrence 
relationships, may be used to solve linear Diophantine equations. This 
method does not require the storage of all the partial results in the Euclidean 
Algorithm. Given 


P-2 = ° 

p.j = 1 

*-2 = 1 

'Ll = 0 

These integer values are the initial conditions of the following respective 
recurrence relationships, where k = 0: 

r k-2 " r k-l h k + r k 

i“H 

1 

M 

u 

V 

X 

fn 

VI 
O 


Pk-l^ + Pk-Z 

(12) 

tl’V + q k-2 

(13) 
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The first expression is the Euclidean Algorithm, from which successive 
values of h^ are determined. In turn, successive values of p^. and are 
determined; the process is terminated at n, where r^ = 0. The final values 
are h^ p^ ^ , and q^ and the linear Diophantine equation has the form 


r ->9 1 

- 2 m - 1 


The quotients P^Al^ are known as successive convergents of the continued 
fraction form of the Euclidean Algorithm (Refs. 3 and 4). 

Returning to example 8, hp through h^ = h^ ^ were computed in 
accordance with (11). Concurrently, (12) and (13) can be applied to deter- 
mine pg and q^, as follows: 

Pq = p jhp + p £ = 1 • 3617 + 0 = 3617 

q Q = q -l h 0 + < l- 2 = ° ' 3617 + 1 = 1 

Pj = PqI^ + p i = 3617 • 1 + 1 = 3618 

q l = q 0 h l + = 1 • 1 + 0 = 1 

P2 = P^2 Pq = 1 + 3617 = 7235 


q 2 = q i h 2 + q 0 = 1 • 1 + 1 = 2 

p^ = P 2? 1 3 + Pj = 7235 • 2 + 3618 = 18,088 

q 3 = q 2 h 3 + q l = 2 ' 2 + 1 = 5 

P 4 = < l 3 h 4 + <1 2 = 5 ' 1+2 = 7 

p 5 = p 4 h 5 + P 3 = 25,323 • 2 + 18,088 = 68,734 

q 5 = q 4 h 5 +q 3 =7-2 + 5 = 19 
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Since n - 1 = 5, the process terminates and (14) can be evaluated. 


r -lP5 ‘ r -2 q 5 = ( " 1)6 (r -2’ r -l ) 

64 • 68, 734 - 231, 525 • 19 = 1 
or 

231, 525 (-19) + 64 • 68, 734 = 1 

The results agree with those in example 8, where 

Yj = -19 = 45 mod 64, and z = 68,734 

Note that only r^ 9^ an ^ 9^. ^ need be stored (current values). 

Values of p^ do not have to be computed, since p^ ^ = z, which is not needed 
in evaluating y^. 

Whenever the values of the m. are large (i. e., exceed several hundred), 
the preceding method is most efficient in terms of iterations (time) and 
storage required by a general-purpose computer. 
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APPENDIX B. AN APL PROGRAM FOR THE CHINESE 
REMAINDER THEOREM 


I. AN APL PROGRAM FOR SOLVING SIMULTANEOUS CONGRUENCES 

OVER PAIRWISE RELATIVELY PRIME MODULI 

APL (A Programming Language) is an interactive programming lan- 
guage created by K. E. Iverson. Complex sequential processes may be con- 
cisely described in APL with a minimum amount of self- training. It is 
particularly suited for testing feasibility of algorithms without extensive 
programming experience. Manipulations on entire arrays of operands can 
be efficiently performed in APL (Ref. 5). 

An APL terminal at the Jet Propulsion Laboratory connects via 
telephone data lines to a time- shared IBM 360 model 50 general-purpose 
computer. The work- space capacity is 48K words. 

The statements comprising an APL program entitled 11 PRIMECOUNT" 
appear in Fig. B-l. Upon the request of the program, the user enters the 
number of counters (i. e., number of congruences to be solved simultaneously), 
each modulus, the count associated with each modulus (the a^ associated with 
the rm in Eq. 6, and the clock frequency. The program determines each y. 
needed in Eq. 7 to determine the number of clock pulses x in order to achieve 
the desired count in each counter simultaneously, the minimum number of 
clock pulses (x reduced modulo M), and the time in milliseconds required to 
generate the minimum number of clock pulses. Time in seconds is computed 
to eight significant digits . 

Figure B-2 is the APL program solution of example 7, with a clock 
frequency of 20 MHz assumed. The boxlike symbols followed by a colon are 
points in time where the user is requested to enter parameters. Note that 
the solution of any number of simultaneous congruences over pairwise rela- 
tively prime moduli may be determined by the APL program providing work- 
space capacity is not exceeded and CPU time is acceptable (i. e., cost). 

II. COMPUTATION OF y VALUES 

Moduli whose values do not exceed several hundred are anticipated in 
future timing and control designs. Therefore, the iterative relations given 
in Appendix I were not used in the PRIMECOUNT APL program to determine 
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the values of y.. Instead, y^, starting with a value of 0, is incremented and 
tested to determine whether or not it satisfies 

M.y. - 1 mod m. 
i 7 1 i 

That is, does 

m.l M.y. - 1 (15) 

r i i 

where rm and M. are given and y. is set at 1? If not, y. is incremented to 2, 
and the test is repeated, etc. The loop may have to be traversed as many as 
m. - 1 times before a y. satisfying (15) is found. For large values of m., 
the Euclidean Algorithm for determining successive h values and the 
recurrence relation for determining successive q values are recommended, 
as shown in Appendix A. For example 7, the APL program performs 
45 iterations to determine y^, whereas six iterations on h^ (h^ through h,.) 
and six on q^ (q^ through q,-) are required to determine the same y^ in 
example 8. 

Note that for a given timing and control system organization, the y. 
are calculated once. Any number of sets of a^ may be entered. State- 
ment [28] in the APL program in Fig. B-l asks whether there are any 
additional counts. The user types a yes or no. A yes causes a branch to 
statement [20], preparing the program to accept a new set of a. (i. e., counts). 
A no terminates the program, as shown in Fig. B-2. 
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V 

IPRIMECOUNTl L]]V 

PRIMECOUNT ;M;Y \R\L -,I -,A;F-,T 


Cl] 

'EE TER THE NUMBER OF COUNTERS' 


C23 

n + n 


[3 ] 

• ENTER EACH MODULUS' 


[4] 

M*-[} 


C 5 ] 

■ +ERR0R x \N*p ,M 


[6] 

y+n p o 


C7] 

R+N pO 


[8] 

L<-x /M 


[9] 

1+0 


CIO] 

I+I + 1 


[11] 

i?[ I >£*/./[!] 


[12] 

J -<“0 


[13] 



[14] 

-*-13xi l*tf[J] | J*R[ J] 


[15] 

m w 


[16] 

-»-10xi N*I 


[17] 

'CORRESPONDING Y VALUES ' ;Y 


[18] 

'ENTER CLOCK FREQUENCY IN MEGAHERTZ' 


[19] 

F-n 


[20] 

BR1:' ENTER COUNT ASSOCIATED WITH EACH MODULUS' 


[21] 

/!<-[] 


[22] 

-+ERR0R* \N*p ,M 


[23] 

X+L | +/L \A*L | Y*R 


[24] 

'CLOCK PULSES REQUIRED X= ' ; X ; 1 + ' ; L ; ' *K FOR K= 0,1,2,.. 

1 

• • 

[25] 

'MINIMUM X= ' ; X 


[26] 

T-*-( Ll00000000x5£’~9 +X tFx1000000 )t100000 


[27] 

'MINIMUM X CORRESPONDS TO '\T\' MILLISECONDS' 


[28] 

'ADDITIONAL COUNTS?' 


[29] 

-+BRI x i 'Y'cH 


[30] 

-*-0 


[31] 

V 

ERROR-.' INVALID INPUT' 



Fig. B-l. An APL program for solving simultaneous congruences over 

pairwise relatively prime moduli (Chinese Remainder Theorem) 
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PRIME COUNT 

ENTER THE NUMBER OF COUNTERS 

□ : 

4 

ENTER EACH MODULUS 

□ : 

64 27 25 343 

CORRESPONDING Y VALUES 45 13 19 19 

ENTER CLOCK FREQUENCY IN MEGAHERTZ 

□ : 

20 

ENTER COUNT ASSOCIATED WITH EACH MODULUS 

□ : 

37 17 2 341 

CLOCK PULSES REQUIRED X- 9 0 3 9077 + 1 4 8 1 7 60 0 *K FOR K- 0 , 1,2 
MINIMUM *=9039077 

MINIMUM X CORRESPONDS TO 451.95385 MILLISECONDS 
ADDITIONAL COUNTS? 

NO 

Fig. B-2. APL solution of example 7 
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